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1.  Introduction 

1.A  Background: 

Each  year  about  1 86,000  women  contract  breast  cancer,  of  whom 
46,000  will  die  of  it.  Breast  cancer  is  the  leading  cause  of  cancer 
deaths  amongst  American  women  aged  35  to  50  and  also  has  the 
highest  cancer  mortality  rate  for  younger  women.  Mammography 
is  most  widely  used  modality  for  early  detection  of  breast  cancer. 

The  detection  of  tiny  carcinomas  and  low  contrast  tissue  is  essen- 
tial  for  mammography  to  reduce  the  mortality  rate  ’  .  It  has  been 
shown  that  about  20%  of  detectable  cancers  are  overlooked  or  mis-  figure  j  Compkte  optic  con. 
diagnosed  on  first  detection.  We  obviously  need  a  major  improve-  structed  bv  confieurinp  fibers 
ment  in  the  quality  of  mammograghic  imaging.  Using  an  innova-  tbrougb  metal  grids. 
tive  new  technology  x-ray  optics  to  improve  system  resolution,  en¬ 
hance  contrast  and  reduce  dosage  will  increase  the  effectiveness  of  this  proven  screening  modality,  with 
a  direct  impact  on  mortality. 

1.B  Benefits  of  digital  mammography 


Conventional  film  screen  mammography  suffers  from  limited  dynamic 
range  and  film  granularity,  which  can  reduce  the  sensitivity  to  detect  micro¬ 
calcifications.  The  primary  theoretical  limitations  of  mammography  are  the 
system  resolution,  which  determines  the  minimum  size  of  the  detectable 
malignancy,  and  the  need  to  expose  the  patient  to  ionizing  radiation.  Using 
innovative  new  technology  —  polycapillary  optics  —  to  improve  system 
resolution  and  reduce  required  dose  will  increase  the  effectiveness  of  this 
proven  screening  modality,  with  a  direct  and  immediate  impact  on  mortality. 
By  matching  polycapillary  optics  to  digital  detectors,  we  can  develop  direct 
digital  mammographic  systems  which  can  avail  themselves  of  the  advan¬ 
tages  of  digital  processing,  including  improved  image  contrast  and  resolu¬ 
tion  at  reduced  radiation  dose.  In  practice,  mammographic  imaging  is  often 
limited  by  quality  assurance  issues,  which  can  also  be  favorably  addressed 
by  digital  processing.  Digital  detection  can  provide  high  dynamic  range 
which,  in  addition  to  improving  contrast,  greatly  increases  the  tolerance  of  the  final  image  to  under  or 
over  exposure.  Digital  images  can  be  enhanced  and  are  amenable  to  computer-aided  diagnosis.  Finally, 
digital  images  can  also  be  transported  quickly  for  skilled  consultations. 


Figure  2  Cross  section 
of  polycapillary  fiber. 
0.5  mm  in  diameter. 
SEM photo. 


1. C  What  capillary  x-ray  optics  can  do 

Capillary  optics,  consisting  of  arrays  of  hollow  glass  tubes,  is  a  relatively  new  technology  for  controlling 
X-ray  beams.  X-rays  incident  on  the  interior  of  the  glass  tubes  at  small  angles  are  guided  down  the  tubes 
by  total  external  reflection.  Arrays  of  curved/tapered  capillaries  can  be  used  to  focus,  collimate  and  filter 
X-ray  radiation3’4,5’6’7.  Such  arrays  are  manufactured  by  stringing  hollow  glass  polycapillaries  through 
metal  grids  (as  in  Figure  1)  or  as  a  monolithic  optic. 
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The  use  of  Kumakhov  capillary  optics  in  place  of  conventional  scatter  reduction  grids  in  a 
mammographic  system  has  significant  potential  to  provide  improved  resolution,  increased  contrast  en¬ 
hancement  and  reduced  dose  in  mammographic  imaging.  The  optics  can  also  be  used  to  mate  the 
radiographic  image  with  a  digital  detector  by  appropriate  choice  of  demagnification  and  separation  to 
discrete  chips.  In  addition,  a  pre-patient  optic  could  be  employed  to  increase  the  available  intensity  in  a 
fan  beam  relative  to  simple  slot  collimating. 

The  fibers  in  this  research  varied  from  0.3  mm  to  4  mm  in  outer  diameter,  with  channel  sizes  (the 
diameters  of  the  holes  within  the  fibers)  ranging  from  4  to  22  pm,  as  shown  in  Table  1  in  section  2.E. 
Open  area  in  the  table  refers  to  the  fraction  of  the  front  face  of  the  fiber  which  is  open  channels,  as 
opposed  to  glass  walls.  No  fiber  can  have  a  transmission  that  is  greater  than  the  open  area. 

X-ray  beams  can  be  bent,  focused  or  collimated  by  a  carefully  curved  polycapillary  array.  X-ray  photon 
energies  are  much  larger  than  the  plasma  frequencies  of  glasses,  which  are  tens  of  electron  volts.  In  this 
regime,  the  real  part  of  the  index  of  refraction  of  glass  can  be  simply  approximated  by 

,  ®  l  (1) 


where  n  is  the  index  of  refraction,  s  is  the  dielectric  constant  of  the  glass,  to  is  the  photon  frequency,  cop 
is  the  plasma  frequency  of  the  material,  and  so  is  the  dielectric  constant  for  vacuum. 

The  plasma  frequency  cop  for  glass  is  small  compared  to  the  photon  frequency,  so  that  n  is  slightly  less 
than  unity.  Therefore,  x  rays  traveling  in  vacuum  or  air  can  be  totally  externally  reflected  from  smooth 
glass  surfaces.  The  maximum  angle  between  the  x-ray  beam  and  glass  surface  for  total  reflection,  0C,  is 
approximately 

0  c  =  mrad,  (2) 

where  E  is  the  photon  energy  in  keV.  The  higher  the  energy  of  photon,  the  smaller  is  the  critical  angle. 
For  this  reason,  profile  errors  that  change  the  angle  of  incidence  are  more  detrimental  at  high  energies. 

1.D  Numerical  simulations 

To  evaluate  the  experimental  performance  of  capillaries  and  optics,  and  to  design  capillary  optics,  it  is 
necessary  to  be  able  to  predict  theoretical  behavior  for  complex  geometries.  Extensive  modeling  pro¬ 
grams  describing  the  propagation  of  x-rays  along  capillaries  and  optics  are  being  developed.  These 
modeling  programs  provide  essential  information  on  the  transmission  efficiency  and  divergence  of  cap¬ 
illary  optics,  and  are  being  used  to  simulate  the  performance  of  the  optics  by  obtaining  the  following 
outputs:  transmission  efficiency,  exit  divergence,  output  uniformity,  and  transmission  as  a  function  of 
entrance  angle. 

The  models  for  these  simulation  programs  are  based  on  Monte  Carlo  simulations  of  simple  geometrical 
optics.  The  simulations  randomly  generate  millions  of  x-ray  photons  from  a  particular  location  and  solid 
angle,  and  then  trace  them  through  the  optics  using  reflection  theory  and  classic  mechanics.  The  neces¬ 
sary  confidence  in  the  simulations  is  being  built  by  comparing  the  results  of  the  modeling  program  with 
experimental  data. 


For  single  fibers,  the  computational  speed  is  greatly  enhanced  by  a  reduction  to  two  dimensions  by  pro¬ 
jecting  the  trajectory  onto  the  local  fiber  cross-section.8  Reflections  are  computed  from  standard  tables.9 
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1.E  Purpose 

The  purpose  of  this  project  is  to  contribute  to  an  ongoing  effort  to 
use  the  new  innovative  technologies  in  X-ray  optics  to  develop  a 
digital  mammographic  imaging  system  with  high  resolution  and 
contrast  by  extending  the  capabilities  of  the  computer  simulation. 

Enhanced  computer  simulations  are  critical  because  manufactur¬ 
ing  large  polycapillary  optics  is  difficult,  expensive  and  time- 
consuming,  and  simulation  is  a  good  design  tool  to  predict  the 
performance  of  the  optics  and  make  it  clear  how  all  the  different 
parameters  work  together  and  impact  the  resultant  transmission, 
resolution  and  contrast. 

1. F  Technical  objectives 

The  specific  aims  of  this  traineeship  are: 

•  to  simulate  and  measure  trial  optic  setups,  to  compare  the 
simulated  results  with  experimental  data  and  to  determine  the 
best  design  for  building  the  best  direct  digital  mammographic 
system; 

•  to  further  the  author’s  knowledge  about  the  mammographic  systems  and  to  promote  a  future  career 
in  this  field. 

2.  Body:  Methods  and  Results 

2. A  Introduction  to  profile  and  surface  defect  corrections 

Real  glass  surfaces  are  not  perfect.  In  addition  to  absorption,  there  is  roughness,  which  has  short  wave¬ 
lengths,  of  the  order  of  microns  or  less,  waviness,  which  has  longer  wavelength,  and  bending,  which  has 
wavelength  longer  than  the  length  of  the  capillary.  Bending  is  described  by  the  radius  of  curvature  of 
the  fiber,  R  in  Table  1 .  Waviness  results  in  the  channel  wall  being  tilted  with  respect  to  the  axis  of  the 
fiber,  and  therefore  changes  the  angle  of  incidence.  It  is  described  in  terms  of  a  distribution  of  tilt  an¬ 
gles.  Surface  roughness  produces  both  diffuse  scattering,  for  which  the  reflected  angle  can  be  either 
larger  or  smaller  than  the  corresponding  incident  angle,  and  transmission  into  the  glass.  The  roughness 
correction10  used  in  previous  simulations11  is  described  in  terms  of  a  roughness  height,  z  and  a  correla¬ 
tion  length,  s. 

2.B  Statement  of  Work 

To  provide  a  framework  for  assessing  the  results  of  the  project,  the  original  statement  of  work  for  the 
full  three-year  project  is  reproduced  below. 

Proposed  Statement  OfWprk 

Task  1  Developing  and  testing  the  simulation  software,  (months  1-18) 

Task  2:  Perform  experiment  and  geometric  simulations  on  pre-patient  collimating  optic 
at  a  wide  range  of  energies  (month  19-23) 

2.1  Measure  transmission  as  a  function  of  x-ray  energy. 
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Figure  3.  The  generation  of  a  ran¬ 
dom  point  (of  a  photon)  in  a  circle. 
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2.2  Measure  transmission  uniformity. 

2.3  Measure  dependence  of  transmission  on  x-ray  entrance  angle. 

2.4  Perform  geometric  simulation  on  transmission  spectra. 

Task  3:  Perform  experiments  on  collimating  optic  at  mammographic  energies. 

(month  24-26) 

3.1  Measure  exit  angle  divergence. 

3.2  Measure  contrast  enhancement. 

Task  4:  Perform  experiments  and  geometric  simulations  on  post-patient  anti  scatter 
optic,  (month  27-30) 

4.1  Measure  central  part  and  whole  optics  as  a  function  of  photon  energy 

4.2  Measure  transmission  uniformity 

4.3  Measure  dependence  of  transmission  on  x-ray  entrance  angle 

4.4  Perform  geometric  simulation  on  transmission  spectra 

Task  5:  Perform  experiments  on  post-  patient  anti  scatter  optic  at  mammographic 
energies  (month  31-33) 

5.1  Measure  transmission  energy  as  a  function  of  photon  energy 

5.2  Measure  transmission  uniformity 

5.3  Measure  transmission  entrance  angle  dependence 

5.4  Measure  scatter  fraction 

5.5  Measure  contrast  enhancement 

5.6  Measure  uniformity  of  magnification  factor 
Task  6.  Design  assessment  (month  34-36) 

In  the  first  year  of  this  project,  we  developed  and  tested  a  more  physical  model  and  compared  the  simu¬ 
lation  results  for  a  variety  of  straight  capillaries  to  measured  data  in  the  energy  range  from  10  to  80  keV. 
Using  this  analysis,  the  number  of  required  parameters  in  this  energy  range,  including  20  keV  for  mam¬ 
mography,  is  reduced  from  four  to  two,  bending  curvature  and  waviness.  It  was  not  found  to  be  neces¬ 
sary  to  include  roughness  effects.  Although  it  took  more  than  a  year  to  develop  this  unexpected  new 
model  and  perform  its  implementation  for  single  fibers  and  lenses,  it  was  worthwhile  because  we  now 
have  a  more  physical  simulation  model  which  increases  the  time  efficiency  for  future  modeling. 

Sections  2C-2E  are  the  comparison  of  the  new  model  with  the  previous  model  and  the  comparison  of 
their  simulation  results  with  the  experimental  data  for  borosilicate  glass  optics.  During  the  second  year 
the  principal  investigation  was  changed  from  Hui  Wang  to  Cari.  Section  2F  contains  measurement  of  a 
post  patient  scatter  rejection  borosilicate  glass  optic.  Section  2G  presents  data  and  simulations  applied 
during  the  second  year  to  lead  glass  optics.  Section  2  H  present  measurements  and  simulations  of  a  long 
focal  length  collimating  optic  applied  during  the  third  year. 

2.  C  Improvements  in  generating  x-ray  photons 

As  shown  in  Figure  3,  to  produce  a  random  point  P  in  a  circle,  we  can  randomly  generate  either  (r,  0)  or 
(x,  y).  Previous  simulations' 1  used  uniformly  generated  9  in  the  range  (0,  2n)  and  r  in  the  range  (0,  R), 
where  R  is  the  radius  of  the  x-ray  source,  capillary  channel  or  fiber  input  area.  This  work  uniformly 
generated  x  and  y  coordinates,  with  the  limitation  x2  +  y2  <  R2,  which  is  closer  to  an  isotropic  distribu¬ 
tion. 
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2.D  Improvements  in  waviness  correction  model 


Waviness  of  a  capillary  surface  is  oscilla¬ 
tions  with  wavelengths  shorter  than  the 
capillary  length  and  longer  than  the 
wavelength  of  the  roughness.  The  wavi¬ 
ness  effect  is  modeled  as  a  random  tilt  of 
the  glass  wall,  which  results  in  a  random 
change  in  the  reflection  angle  on  each 
bounce.  If  the  random  tilt  angle  is  (|>  and 
the  incident  angle  is  0,  then  the  reflection 
angle  with  respect  to  the  nominal  surface 
will  be  0+2(|>,  as  shown  in  Figure  4.  Pre¬ 
vious  simulations  used  a  uniform  distri¬ 
bution  model;  this  work  uses  a  tilt- 
corrected  normal  distribution  to  generate  <|>. 

2.D.1  Uniform  distribution  model 

Previous  simulations,  referred  here  as  Model 
Ml ,  assumed  a  uniform  distribution  of  2<J> 
from  -A0max  to  AQmax  if  0  >  A0max,  or  from  -0 
to  A0max  if  0  <  A0max-  The  maximum  random 
tilt  angle  A0max  is  the  adjustable  parameter 
which  describes  the  amount  of  waviness  of  the 
capillary. 

2.D.2  Tilt-corrected  normal  distribution 
model 

2.D.2.1  Normal  distribution 


Figure  4.  Scheme  of  an  x-ray  beam  interactions  with  a 
randomly  lilted  surface. 
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/ 
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Figure  5.  Three  surfaces,  OA],  OA2  and  OA3,  with 
different  tilt  angles,  <f>],  fc  and  fo,  respectively,  from 
the  nominal  surface  OO  ’.  Here,  -6  <fo<0, 


OAf 


'COS  (j)  1  =  OA2*COS  (/>  2=  OAj*COS  (f>  3  =  OO’. 


For  this  work,  named  as  Model  M2,  it  is  as¬ 
sumed  that  these  tilt  angles,  (J),  are  normally 
distributed  in  the  range  (-ti/2,  7i/2).  with  the 

mean  value  equal  to  zero.  For  high  quality  optics,  the  standard  deviation  of  this  normal  distribution,  <7,  is 
much  smaller  than  the  critical  angle,  0C.  The  probability  distribution  of  tilt  angles,  <|>,  is 


GW)  = 


1  -V_ 

i  o  _  2 


a  -Jljx 


2  <7 


(3) 


2.D.2.2  X-ray  impact  tilt  correction 

In  this  work,  Model  M2,  consideration  was  taken  of  the  fact  that  the  surface  tilt  angle  will  affect  the 
probability  of  x-ray  impact  on  that  surface.  Taken  to  extremes,  a  surface  region  perpendicular  to  the 
beam  is  much  more  likely  to  intercept  the  beam  than  a  surface  region  parallel  to  the  beam.  Figure  5  dis¬ 
plays  three  surfaces,  OAi,  OA2  and  OA3,  with  different  tilt  angles,  <|>i,  <J>2  and  4)3,  respectively,  from  the 
nominal  surface  OO’.  The  projections  onto  the  nominal  surface  for  the  three  surfaces  are  equal, 
OAi*cos<|>i  =  OA2*cos<j>2  =  OA3*cos<))3  =  OO’.  The  probability  of  incidence  of  a  parallel  x-ray  beam  with 
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incident  angle  0  (with  respect  to  the  nominal  surface  00’)  hitting  these  tilted  surfaces  is  given  by  their 
corresponding  perpendicular  length  AjBj  (for  j  =  1,  2  or  3) 

P.  oc A,B.  =  0At  ■  sin(# +4>j)=  ■  sin((9  +  <j>; ), 

COStpj 


We  call  this  the  tilt-corrected  probability  distribution.  The  complete  description  of  this  distribution  is 


Fsin(g  + 
cos(^) 

0, 


(5) 


where  0  is  the  incident  angle,  <|>  is  the  tilt  angle  and  Fisa  normalization  constant. 


2.D.2.3  Tilt-corrected  normal  distribution 


Combining  the  normal  distribution  G(<j>)  with  the  tilt  correction  H(<|))  gives  <j>  for  a  certain  incident  angle 
0  as  follows: 


RP  sin(<9+^ 
cos^ 

0, 


- e«f>< | 


where  0  is  the  incident  angle,  §  is  the  tilt  angle  and  K  is  a  normalization  constant. 
Noting  that 

sin(9+^)  _sjnfl+cosfl.tan^ 
cos^ 


(6) 


(7) 


and  that  both  the  incident  angle,  0,  and  the  tilt  angle,  <j),  are  very  much  less  than  40  mrad,  we  use  the  ap¬ 
proximation 


-r 

Ce 2al  *(0 +  </>), 

0, 


-9  <<j>  < 


n 


(8) 


it 


--<<f><-9 

2 


where  the  normalization  constant ,  C,  is 


C 


1 

-s/2^r  9a  +  2  a 2 


(9) 


to  speed  the  calculation.  P((|))  is  called  the  tilt-corrected  normal  distribution.  In  model  M2,  a  is  the  vari¬ 
able  that  describes  the  amount  of  waviness.  Two  examples  of  J((j>)  and  P(4>)  are  shown  in  Figure  6  and 
Figure  7. 
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Fiber  Description 

Model  Ml 

Model  M2 

Fiber 

# 

Type 

Outer 

Diame¬ 

ter 

mm 

Channel 

Size 

pm 

Open 

Area 

Length 

mm 

z 

nm 

s 

pm 

R 

m 

AOmax 

mrad 

R 

m 

CT 

mrad 

A 

1 

0.5 

12 

65% 

105 

0.7 

6 

105 

0.4 

60 

0.225 

C 

3 

0.75 

22 

50% 

136 

0.5 

6 

125 

0.35 

225 

0.2 

D 

4 

4 

12 

55% 

130 

0.8 

6 

110 

0.285 

139 

0.125 

E 

5 

0.3 

4-5 

55% 

105 

0.7 

6 

28 

0.2 

31 

0.09 

F 

4 

4 

12 

55% 

130 

0.8 

6 

90 

0.45 

90 

0.18 

Table  1  Parameters  for  best-fit  simulations.  R  is  the  bending  radius.  For  Ml,  AQmaxis  the  amount 
of  waviness,  z  is  the  roughness  height  and  s  is  the  roughness  correlation  length.  For  M2,  cris  the 
standard  deviation  of  the  waviness. 


Figure  6.  The  probability  distribution,  J(</>),  and  Figure  7.  J(</>)  and  P(<j>)  versus  tilt  an- 

its  approximation,  P  (</>),  versus  tilt  angle  (j>  when  gle  (j)  when  incident  angle  6  =  0.009  rad,  stan- 
incident  angle  6  =  0.0001  rad,  standard  deviation  dard  deviation  a  =  0.0004  rad.  <f>  is  in  rad. 
u  =  0.0004  rad.  (j)  is  in  rad. 
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X-ray 


Pinhole 


Detector 


2.E  Results 

_ ^  Transmission  of  a  number  of  different  fibers 

I®®  I  Fiber  ' - ^  measured  previously11  in  the  energy  region 

/  A  ^1  . —  10-80  keV  are  shown  in  figures  9-19.  The  de- 

scriptions  of  these  polycapillary  fibers  and 

- ^  /  \  ' - ^  their  simulation  parameters  are  shown  in 

p.  *  Table  1 .  The  experimental  setup  for  straight 

X-ray  1  0  e  !  fibers  measurement  is  shown  in  figure  8.  A 

Source  Lead  Detector  microfocus  source  with  a  tungsten  target  and 

shielding  source  spot  size  of  50  pm  was  used  in  the  ex¬ 

periment.  The  source  was  mounted  on  two 
Figure  8.  Experimental  setup  for  single  fiber  stages,  which  could  be  moved  in  two  direc- 
measurement  tion  transverse  to  the  beam,  and  also  enclosed 

in  a  thick  lead  box  to  reduce  scattered  x-rays. 
The  distance  between  the  source  and  the  fiber  was  about  100  cm  to  ensure  incident  angles  less  than  criti¬ 
cal  angle  up  to  energy  80  keV.  A  200  pm  pinhole  was  attached  to  thick  lead  shielding.  The  size  of  the 
pinhole  is  smaller  than  the  diameter  of  the  fiber  which  is  about  0.5  mm  to  prevent  scattering  and  leak¬ 
age,  to  reduce  the  detector  dead  time  and  to  keep  the  distance  between  the  source  and  the  fiber  as  small 
as  possible.  The  fiber  was  put  on  a  finely  machined  groove  in  an  aluminium  plate  mounted  on  two 
stages,  which  could  be  moved  in  two  direction  transverse  to  the  beam. 


Lead 

shielding 


Figure  8.  Experimental  setup  for  single  fiber 
measurement 


*  data  ^  V'0 

— Ml:  R=80  ■" 

-  -o-  ■  WE  R=80  \'| 

-#-•  WH:  R=100  \D 

-o-M2:R=100  \ 


o  u  □- 
A  '  K  ‘S' 


•c  X  data  * 

I35-  - T- •  * 

p  M2o=0.15 

30.  -  o  IVEa=0.25 

a  IV^aO.3  5 


*  data  . 
— M2:  W=0.3, 1000  photons  “ 

-  -o-  -  M2:  W=0.3, 10000  photons 

-  a  -  M2:  W=0.4,1000  photons 
--O-  -  W12:  W=0.4, 1 0000  photons  * 


20  40  60  80 

Photon  energy(Ke\/) 


20  40  60  80 

Fhotonenag^l^ 


60  80 


Photon  energy(KeV) 


Figure  9  Transmission  spectra  Figure  10  Simulations  of  Figure  11  Simulations  of  trans- 

for  fiber  C  simulated  with  transmission  spectra  for  fiber  mission  spectra  with  waviness 

different  bending  curvature  C  with  only  waviness  com-  correction  using  different  numbers 
alone  using  models  Ml  and  pared  with  the  experimental  of  photons,  compared  with  the  ex- 

M2,  compared  with  expert-  data.  This  figure  shows  the  perimental  data  of fiber  C.  The 

mental  data.  This  figure  shows  effects  of  waviness.  It  also  figure  shows  that  the  simulation 

that  the  simulations  with  shows  simulations  using  curves  using  10000 photons  are 

bending  alone  can  not  fit  the  waviness  alone  do  not  fit  the  much  smoother  than  those  using 

data,  and  indicates  that  data.  The  simulations  do  not  1000  photons  because  of  the  re¬ 
bending  is  not  the  only  factor  include  the  roughness  or  duclion  in  statistical  fluctuation, 

to  cause  the  transmission  bending  correction, 

drop. 
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The  experimental  data  and  simulation  results  from  the 
two  simulations  for  a  type  3  polycapillary,  Fiber  C,  are 
compared  first  to  demonstrate  the  fitting  processes,  to 
interpret  the  waviness  and  bending  effects,  and  to  un¬ 
derstand  the  final  results. 

2.E.1  Bending  effects 

A  slight  bending  can  dramatically  reduce  the  transmis¬ 
sion  of  high-energy  photons  because  of  the  small  criti¬ 
cal  angle  at  high  energies.  Figure  9  shows  a  compari¬ 
son  between  experimental  data  and  simulations  with 
different  bend  radii.  The  difference  in  the  bending  ef¬ 
fect  between  Ml  and  M2  is  due  to  the  different  algo¬ 
rithms  used  in  generating  x-ray  photons  as  discussed  in  Figure  12  Simulated  transmission  spectra 
section  0.  Ml  assumes  more  photons  hit  the  center  of  using  modd  M2  with  fixed  waviness  (a=0.2 
the  fiber  as  well  as  the  center  of  each  channel.  M2  gen-  mmd)  and  different  bending,  compared  with 
erates  photons  uniformly  across  the  channel,  fiber  and  the  experimental  data.  The  simulation  curves 
source.  shown  are  just  a  few  representative  cases  in 

For  both  models,  a  bend  of  R  =  80  m  reduces  the  ^e  bisection  process.  For  M2,  the  best-fit 

transmission  too  much  at  high  energies.  Since  any  parameters  are  <7=  0.2  mrad  and  R  =  225 
waviness  or  roughness  correction,  like  the  bending  cor-  m;  for  Ml,  they  are  z  =  0.5  nm,  s  =  6  pm, 
rection,  will  reduce  the  transmission,  R  must  be  larger  A  9max  =  0.35  mrad  and  R  =  125  m. 
than  80  meters.  We  set  Rmj„  =  80  m,  which  is  useful  in  determining  the  best-fit  parameters  for  the  model 


Photon  energy(KeV) 


Source  position(mm) 

Figure  13.  Simulations  of  source  scan  curve  us¬ 
ing  model  M2  with  different  waviness  ( o)  com¬ 
pared  with  the  experimental  data.  The  larger  the 
waviness,  the  narrower  the  simulated  scan  curve. 
The  simulation  curves  with  the  best  fitting  pa¬ 
rameters  of  Ml  and  M2  obtained  in  Figure  12 
are  also  shown  in  this  figure  for  reference. 


-4-3-2-10  1  2  3  4 

Source  position(rmri) 


Figure  14  Simulations  of  source  scan  curves 
using  model  M2  -with  no  waviness  and  differ¬ 
ent  bending  compared  with  the  experimental 
data  at  10  keV.  The  simulated  scan  curve 
does  not  become  narrower  even  when  the 
bending  radius  R  =  20  m.  The  word  “ideal" 
refers  to  a  simulation  with  R-oo  and  no 
waviness. 
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M2.  Clearly,  bending  alone  cannot  explain  the  experi¬ 
mental  data,  especially  at  low  energies. 

2.E.2  Waviness  effects 

To  explain  the  data  at  energies  from  20  to  60  keV,  an¬ 
other  correction  is  needed.  Roughness  does  not  have  a 
major  effect  in  this  energy  range11.  Waviness  is  ex¬ 
pected  to  be  the  major  factor  in  this  energy  range. 

Figure  10  shows  that  the  transmission  becomes  lower 
for  all  the  energies  and  the  high-energy  transmission 
drops  faster  with  increasing  a  gets  larger.  It  also  shows 
that  waviness  alone,  like  bending  alone,  can  not  ex¬ 
plain  the  experimental  data.  In  the  simulation,  a  =  0.3 
mrad  and  0.4  mrad  are  definitely  over-correcting  at  low 
energies,  so  we  can  set  0.3  mrad  as  the  maximum  of  a 
(CTmax)  for  Fiber  C.  The  waviness  parameter  A0  for 
model  Ml  is  larger  than  the  parameter  a  for  model  2 
because  the  tilt  correction  in  model  M2  makes  it  more 
probable  that  the  x-ray  will  intersect  the  high  angle  sur¬ 
faces. 

The  fluctuation  of  the  simulation  curve  in  Figure  1 0 
when  cr=0.25  mrad  is  mainly  due  to  the  relatively  small 


Source  position(mm) 

Figure  15  Simulations  of  source  scan  curves 
using  model  M2  with  fixed  waviness  (<j=  0.2 
mrad)  and  different  bending  compared  with 
the  experimental  data.  The  simulated  scan 
curve  does  not  apparently  change  even  when 
the  bending  radius  R  =  20  m.  With  and  without 


the  bending,  simulations  using  cr  =0.2  mrad fit 


the  experimental  data  quite  well. 


number  of  photons,  1000,  used  in  the  simulation.  Figure  1 1  shows  the  difference  when  the  number  of 


photons  changes  from  1000  to  10000. 


Using  a  bisection  algorithm  with  Rmjn  and  crmax,  R  and  a  were  applied  together  and  compared  to  the  data 
to  get  the  best  fit  parameters,  as  shown  in  Figure  12.  The  best  fit  parameters  for  fiber  C  using  model  M2 
are  R  =  225  m  and  a  =  0.2  mrad. 

2.E.3  Verification  using  the  source  scan  curves 

Figure  13  shows  the  source  scan  curve,  that  is  transmission  as  a  function  of  lateral  source  displacement, 
and  simulation  results,  for  increasing  waviness  for  Fiber  C  at  10  keV.  The  simulated  source  scan  slowly 
becomes  narrower  when  cr  gets  larger,  however  that  change  is  too  small  to  be  used  to  determine  a. 
Figure  14  displays  simulations  for  Fiber  C  at  10  keV  with  different  R  and  no  waviness.  Figure  15  shows 
the  effect  of  different  R  with  a  fixed  a.  From  these  two  figures,  it  is  apparent  that  changing  R  from  20 
meters  to  infinity  does  not  have  any  visible  effect  on  the  fitting  curve.  Therefore,  we  can  not  use  the 
simulations  of  the  source  scan  curves  to  determine  either  a  or  R,  but  they  can  be  used  to  verify  a. 

For  Fiber  C,  as  shown  for  10  keV  in  Figure  13,  at  higher  energies  in  Figure  16,  and  as  a  function  of  en¬ 
ergy  in  Figure  12,  the  simulation  curve  of  a=0.2  mrad  fits  the  experiment  data  very  well  with  one  fewer 
fitting  parameter  than  the  best  fit  of  Ml . 
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20KeV  40KeV:  59KeV:  80KeV: 

*  data  *  data  *  data  *  data 

--■--Ml:  Best  fit  M 1 :  Best  fit  — Ml:  Best  fit  — ■ —  M 1 :  Best  fit 

■  -  O- ■  ■  M2  a=0  2  --o-- M2: c=0.2  --0-- M2:o=0.2  - -O- -  M2:ct=0.2 

-A-- M2:cr=0.2,R=225  --A  ••  M2:o=0.2,R=225  *A  •  ■  M2:o=0.2,R=225  M2:o=0.2,R=225 


-2  0  2  -2  0  2  -2  0  2  -2  0  2 


Source  position(mm)  Source  position(mm)  Source  position(mm)  Source  position(mm) 


Figure  16  Simulated  source  scan  curves  using  Ml  and  M2  compared  with  experimental  data  at  four 
different  photon  energies.  They  all  fit  the  data  quite  well.  For  M2,  this  figure  also  verifies  that  the 
bending  has  no  visible  effect  on  the  source  scan.  The  best-fit  parameters  for  Ml  are  R  =  125  m,  A0max 
=  0.35  mrad,  s  =  6  pm  and z  =  0.5  nm. 

2.E.4  Best-fits  of  all  the  fibers  and  analysis  of  the  results 


Using  the  algorithm  described  above,  the 
best-fit  parameters  of  M2  for  Fiber  A,  D 
and  E  were  obtained.  The  best-fit  parame¬ 
ters  for  both  Ml  and  M2  are  shown  in  Ta¬ 
ble  1.  The  best-fit  curves  of  Ml  and  M2 
along  with  the  experimental  data  are  shown 
in  Figures  17-19.  The  fibers  are  described 
in  Table  1. 

In  Figure  18,  the  transmission  for  fiber  A 
shows  a  rapid  drop  for  energies  above  30 
keV.  Although  fiber  D  has  lower  fractional 
open  area  than  fiber  A,  its  transmission  ex¬ 
ceeds  that  of  fiber  A  at  energies  above  30 
keV.  This  is  because  fiber  A  is  thin  (0.5 
mm  in  outer  diameter)  and  flexible,  there¬ 
fore  difficult  to  keep  straight  in  the  meas¬ 
urement  apparatus.  Model  M2,  which  can 
vary  only  bending  and  waviness,  requires  a 
much  sharper  bend  for  fiber  A  than  for  fi- 


Photon  energy  (KeV) 

Figure  1 7.  Simulations  of  transmission  spectra  of fiber 
C,  fiber  D  and  fiber  E  with  their  best-fit  parameters 
compared  with  the  experimental  data.  All  parameters 
related  to  this  figure  are  listed  in  Table  1. 
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ber  D.  This  bending  is  more  significant  at  high  energies,  where  the  critical  angles  are  smaller. 

Transmission  curves  for  fiber  C  and  D  are  similar  in  shape  and  are  nearly  flat  up  to  60  keV,  as  shown  in 
figure  16.  However,  the  smaller  channel  size  for  fiber  D  results  in  the  simulation  being  less  sensitive  to 
bending.  Therefore  a  larger  bending  curvature  (smaller  radius)  is  required  for  fiber  D.  The  waviness 
correction  for  fiber  D  must  therefore  be  smaller  to  keep  the  transmission  about  the  same  as  for  fiber  C. 

Fiber  E  is  the  thinnest  fiber.  Table  1  shows  that  it  has  the  largest  bending  correction  because  of  its  flexi¬ 
bility.  Its  transmission  curve  is  flat  up  to  40  keV.  This  is  due  to  its  very  small  channel  size,  as  discussed 
in  the  comparison  of  fiber  C  and  fiber  D.  However,  if  the  channel  size  is  too  small,  as  in  fiber  E,  it  also 
results  in  more  reflections  being  needed  for  a  photon  to  traverse  the  fiber  and  may  have  introduced  other 
defects  such  as  blocked  channels.  This  is  why  the  transmission  is  only  40%  for  the  energies  below  40 
keV  although  the  open  area  is  around  55%. 


The  high  transmission  and  the  simulation  results  show  that  the  quality  of  the  capillary  fibers  is  quite 
good.  The  bending  radius  is  above  130  meters 


for  type  3  and  type  4  capillaries.  It  is  hopeful 
that  we  can  further  improve  the  high-energy 
transmission  performance  of  polycapillary  fi¬ 
bers  by  decreasing  to  an  optimum  channel  size, 
and  making  them  more  rigid. 

Because  Model  M2  is  physically  plausible, 
gives  reasonable  results,  and  matches  the  ex¬ 
perimental  data  very  well,  as  shown  in  figures 
17  and  18,  we  can  use  it  to  simulate  other  poly- 


0  20  40  60  80 


Photon  energy  (KeV) 

Figure  18.  Simulations  of  transmission  spectra  of 
fiber  A  and  fiber  D  with  their  best-fit  parameters 
compared  with  the  experimental  data.  All  parame¬ 
ters  related  to  this  figure  are  listed  in  Table  1. 

capillary  optics.  Compared  to  Ml,  which  uses  four 
parameters,  M2  uses  only  two  parameters,  which 
makes  simulating  the  experimental  results  more 
straightforward,  easier  and  timesaving.  In  addition, 


Photon  energy  (keV) 


roughness  does  not  have  much  effect  in  the  energy 


range  10-80  keV,  so  M2  gives  us  more  confidence 


Figurel9  Simulations  of  transmission  spectra  w’^1  *'ewer  Parameters- 


of  fiber  D  (unexposed)  and  fiber  F  (exposed)  2.E.5  An  analysis  of  a  radiation-damaged  fiber 

with  their  best-fit  parameters  compared  with 

the  experimental  data.  Layer  =  35  means  that  The  measured  transmission  spectra  and  source 
there  is  a  35  micron  glass  layer  added.  All  pa-  scans  provide  sensitive  tools  for  analyzing  capillary 
rameters  related  to  this  figure  are  listed  in  Ta-  ejual  i  t  .  This  method  was  used  to  analyze  the  effect  of 
ble  1  except  the  glass  layer.  radiation  damage  on  fiber  F  that  has  the  same  type  and 

same  length  as  fiber  D.  Radiation  damage  has  been 
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observed  after  a  very  large  dose  of  x-ray  exposure.12  It  was  found  that  prolonged  exposure  to  synchro¬ 
tron  white  beam  radiation  could  cause  a  measurable  bend  for  a  thin  fiber.  This  bending  can  be  reduced 
by  annealing  and  by  holding  the  fiber  rigid  during  exposure.  Fiber  F  had  1.8  MJ/cm2  exposure  before  the 
current  measurements.  There  is  no  visually  observable  bending  or  other  visible  damage. 

The  transmission  spectrum  and  its  simulation  for  fiber  F  are  compared  with  an  unexposed  capillary 
of  the  same  type,  fiber  D,  in  Figure  17.  The  measured  data  shows  a  significant  transmission  drop  above 
30  keV  and  below  10  keV.  The  simulation  indicates  that  the  exposed  capillary  has  more  bending  and 
waviness.  The  bending  radius  was  changed  from  139  m  to  90  m  and  the  waviness  a  changed  from  0.125 
mrad  to  0.1 8  mrad.  Although  these  profile  changes  have  a  significant  effect  the  transmission,  they  are 
still  very  small  and  not  visible.  The  transmission  drop  below  20  keV  can  be  explained  by  adding  35  //m 
x-ray-absorbing  glass  layer  as  shown  in  Figure  19. 13 

2.F  Measurement  of  a  prototype  post-patient  monolithic  optic 

Polycapillary  optics  can  be  used  for  post  patient  scatter  rejection.  Off-axis  rays  are  not  transmitted 
down  the  optic,  as  shown  in  figure  15.  Monolithic  optics  are  shaped  from  a  single  large  fiber.  The  ba- 
sicgeomtric  of  this  post-patient  monolithic  optic  as  a  cone  with  length  of  the  optic  =129  mm,  output  di¬ 
ameter  8.25  mm,  input  diameter=6.1  mm,  channel  size  =  25pm  and  focal  distance  =  500  mm. 


Figure  20.  Measured  acceptance  angle  as  func¬ 
tion  of  source-optic  distance  for  the  monolithics 
optic  at  20  keV 


Figure  21.  Measured  transmission  vs  source-optic 
distance  at  20  keV. 


The  x-ray  generator  used  in  the  experiment  is  the  same  as  for  the  previous  section.  Variation  of  ac¬ 
ceptance  angle  of  the  optics  as  a  function  of  source  optic  distance  was  used  primarily  to  determine  both 
focal  distance  of  prototype  monolithic  optic  and  the  input  depth  field.  The  acceptance  angle  is  the  full 
width  half  maximum  of  the  source  scan  divided  by  the  source  to  optic  distance.  The  full  width  half 
maximum  will  describe  the  maximum  diameter  of  a  source  over  which  optic  will  collect  efficiently  and 
the  intensity  transmitted  through  the  optic  follows  a  gaussian  curvature  as  shown  in  figure  22.  The  ac¬ 
ceptance  angle  is  a  minimum  at  the  focal  distance,  where  the  transmission  is  maximized.  Focal  distance 


Wang,  Hui  and  Cari 


17 


Study  of  Capillary  X-ray  Optics  for  Digital  Mammography 


of  this  optic  was  determined  by  moving  the  source  in  x  and  y  direction  to  maximize  the  transmission  at 
certain  position  of  the  optic.  The  highest  transmission  as  the  input  focus  distance,  and  the  full  width  at  3A 
of  the  maximum  is  defined  as  the  depth  of  field. 


Figure  22.  Measured  scan  source  at  the  focal 
distance. 


Energy  (keV) 


Figure  23.  Transmission  versus  energy  at 
the  focal  distance 


Figure  21  shows  that  the  focal  length  and  the  depth  field  of  this  optic  are  500  mm  and  92  mm,  re¬ 
spectively.  This  focal  length  is  in  good  agreement  with  the  focal  length  determined  by  using  the  trans¬ 
mission  as  function  of  source  to  optic  distance,  as  shown  in  figure  21.  The  high  angle  transmission  of 
this  optics  close  to  zero,  as  can  be  seen  in  figure  22  when  the  source  is  moved  off-axis.  This  will  give 
good  scatter  rejection.  Finally,  the  transmission  efficiency  of  this  optic  was  measured  at  the  focal  point, 
as  in  figure  23.  However  the  transmission  of  the  long  optic  is  poor. 

2.F  Measurement  and  analysis  of  leaded  glass  fibers 

Leaded  glass  polycapillary  optics  could  be  a  candidate  for  scatter  rejection  because  of  lead’s  superior 
absorption.  Lead  glass  optics  could  be  shorter  and  therefore  easier  to  manufacture  with  better  transmis¬ 
sion.  Lead  glass  optics  may  be  the  fast  way  to  be  bring  capillary  optics  to  clinical  use.  In  this  section,  we 
will  analyze  the  experimental  results  of  a  number  of  leaded  glass  single  capillaries  using  model  M2. 
Transmission  of  a  leaded  glass  polycapillary  in  the  9-80  keV  energy  range  is  shown  in  figure  24.  The 
description  of  these  capillaries  and  the  simulation  parameters  are  listed  in  table  II.  The  characteristics  of 
the  transmission  spectra  of  a  leaded  glass  fiber  are  similar  to  the  radiation-damaged  borosilicate  fiber, 
i.e,  the  transmission  drops  off  at  both  lower  energies  and  higher  energies,  as  shown  in  figure  24.  In  ad¬ 
dition  to  the  waviness  and  bending  simulation  parameters,  the  simulation  has  to  include  a  leaded  glass 
layer  of  33  pm  thick  and  roughness  1 .7  nm  to  fit  the  experimental  data. 
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Fiber 

Simulation 
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ber 

type 

Outer 

Di¬ 

ameter 

mm 

Chan¬ 

nel 

size 

pm 

Open 

area 

% 

Lengt 

h 

mm 

Bend¬ 
ing  R 
m 

Wavi¬ 
ness  w 

mrad 

Laye 

rt 

pm 

Rough¬ 
ness  z 

nm 

A 

0.53 

12 

60 

60 

40 

0.26 

45 

1.8 

B30 

0.51 

12 

60 

30 

30 

0.15 

18 

1.7 

B60 

0.51 

12 

60 

60 

56 

0.15 

33 

1.7 

C 

0.51 

11 

50 

60 

27 

0.19 

33 

1.8 

D 

0.53 

11 

50 

60 

28 

0.20 

43 

1.8 

E 

0.53 

11 

40 

60 

45 

0.28 

47 

2.1 

Table  II.  Fiber  description  and  simulation  parameters 


Energy  (keV) 

Figure  24.  Simulated  transmission  with  and 
without  leaded  glass  filter  for  type  B60  fiber  are 
compared  with  the  experimental  data.  Best  fit  is 
simulated  transmission  with  R=56  m,  w=0.15 
mrad,  z=l.  7  nm,  and  t=33pm. 


%  exp:  9  keV 


Source  Position  (mm) 

Figure  25.  Simulation  of  source  scan  curves  with 
z=l.  7  nm ,  R=56  m,  and  w=0. 15  mrad,  and  with 
roughness  z=l.  7  nm  alone,  compared  with  the  ex¬ 
perimental  data.  Since  the  actual  transmission  at 
this  energy  is  quite  low,  the  source  scan  curve 
might  be  somewhat  distorted,  and  thus  the  simula¬ 
tion  does  not  fit  the  data  well. 
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Figure  26.  Simulated  scan  curves  compared  with  the 
experimental  data  at  three  photon  energies. 


Energy  (keV) 


Figure  27. The  simulation  of  scan  curves  with 
best  fit  parameters  at  50  and  70  keV,  com¬ 
pared  with  the  experimental  data. 


The  best  fit  simulation  compared  with  the 
experimental  data  for  source  scan  curveare 
shown  in  figure  25-26.  At  9  keV  the  simulated 
for  source  scan  does  not  fit,  as  shown  in  figure 
25.  This  may  be  partly  because  the  actual  value 
of  the  transmission  at  9  keV,  as  shown  in  figure 
24,  is  quite  low.  But  the  simulation  at  high 
energies  for  source  scan  curve  were  fit  the  data 
fairly  well. 

Finally,  the  measured  transmission  spectra  for 
a  number  of  leaded  glass  fibers  II  are  shown  in 
figure  28. 

The  models  fit  the  transmission  spectra  well, 
with  the  parameters  listed  in  table  II. 


Figure  28  Simulation  of  transmission  spectra  of ' fi¬ 
ber  A,  fiber  B60  (B  with  length  of  60  mm)  and  fiber 
B30  (B  with  length  of  30  mm)  with  their  best  fit  pa¬ 
rameters  listed  in  Table  II,  compared  with  the  ex¬ 
perimental  data. 
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2.H.  Measurements  and  Analysis  of  a  Long  Focal  Length  Collimating  Lens 


-3-2-10  1  2  3 

Source  Position  (mm) 


Figure  29.  Measured  source  scan  curve 
from  20  -80  keV.  The  figure  shows  that 
there  was  no  x-rays  leakage. 

area  digital  x-ray  detector. 


Another  use  of  poly  capillary  optics  is  beam  shaping, 
collimation,  image  demagnification  were  measure  using 
a  multifiber  borosilicate  optic.  Each  single  fiber  is  about 
0.51  mm  in  diameter,  and  contains  hundreds  or 
thousands  capillary  tubes,  10  pm  in  diameter.  The 
multifiber  polycapillary  optic  was  constructed  by 
stringing  3039  flexible  single  fibers  through  metal 
holding  grids.  The  fibers  at  one  end  of  the  lens  are 
parallel,  and  at  other  end  are  pointed  toward  a  common 
focal  point,  2  m  away.  This  means  that  this  lens  can  be 
used  to  collimate  or  focus  x  rays.  In  collimating  mode, 
the  lens  is  used  as  a  collimator  to  collect  x-rays  from  a 
point  source  at  the  focus  and  produce  a  nearly  parallel 
beam.  In  focusing  mode,  the  lens  is  used  as  a  collector 
for  parallel  a  beam  input  and  produce  a  focused  beam. 
The  focused  beam  can  be  used  to  demagnify  an  image. 
This  is  of  potential  benefit  for  direct  imaging  with  small 


The  focal  length,  transmission  efficiency,  uniformity,  and  exit  divergence  were  measured  in  collimating 
mode.  Scatter  rejection,  contrast  enhancement,  and  demagnification  were  demonstrated  in  focusing 
mode. 


2.H.1  Single  Fiber  Measurements  and  Analysis 
2.H.l.a  Straight  Fiber 

To  ensure  that  source,  pinhole,  fiber  and  detector  were  aligned  well,  the  transmission  as  a  function  of 
source  location  was  measured,  as  shown  in  figure  29  for  a  fiber  with  length  of  140  mm.  The  figure 
shows  that  there  was  not  any  x-ray  leakage  detected.  Then  the  transmission  of  the  fiber  type  used  to  con¬ 
struct  a  collimating  lens  were  measured  in  the  10-80  keV  energy  range.  Several  fibers  of  different 
lengths  were  measured,  three  fibers  from  which  are  presented  in  figure  30.  The  fractional  open  area  of 
the  fiber  is  70  ±  3  %. 
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2.H.l.a.l  Simulation  Analysis 


Figure  30.  Source  scan  with  length  of 
Jibes  is  140  mm.  The  lines  are  simulated 
data  with  waviness,  w=0.12  mrad, 
roughness,  z=0.65  nm,  and  bending, 
R-210  m 


Energy(keV) 


Figure  31.  Transmission  versus  energy 
for  single  fibers  with  three  different 
lengths.  The  lines  are  simulated  data 
with  waviness  w=  0.12  mrad,  and 
roughness  z=  0. 65  nm,  but  with  differ¬ 
ent  bending  curvature,  R] 40=2 10  m, 

R  125=105  m,  R6o=48  m. 


Figure  32.  Simulated  transmission  com¬ 
pared  with  the  experimental  data  for  a 
defected  a  fiber.  Best  fit  is  simulation 
with  waviness  w=0. 12  mrad,  roughness 
z=0.65  nm  bending  radius  R=54  m,  and 
glass  layer  thickness,  t=75  pm. 

energies.  All  three  lengths  of  fibers  shown 


The  experimental  data  then  was  analyzed  using  the 
geometrical  optic  simulation  model  M2  to  understand  the 
performance  of  the  fibers.  The  best  fit  simulation 
parameters  were  determined  in  the  manner  of  section  2  E. 
Using  waviness  from  the  source  scan  curve  at  30  keV, 
roughness  from  source  scan  curve  at  10  keV  and  adding  a 
bending  radius  of  210  m,  a  fit  was  achieved  across  the 
whole  spectrum.  The  simulated  source  scan  curves  with 
these  three  best  fit  simulation  parameters  for  a  140  mm 
long  fiber  compared  with  the  experimental  data  is  shown 
in  figure  30.  Then  the  fitting  processes  were  applied  to 
two  fibers  with  length  of  125  mm  and  60  mm.  The 
simulated  transmission  spectra  for  these  three  fibers 
compared  with  experimental  data  are  also  shown  in  figure 
31.  Each  simulation  fits  the  experimental  data  well.  Since 
the  fibers  are  thin,  a  shorter  fiber  bends  more  and 
therefore  transmits  somewhat  more  poorly  at  high 
in  figure  31  have  good  transmission  up  to  energy  50  keV. 


2.H.l.a.2  Defected  Fiber 


One  of  the  measured  fibers  contained  defects  which  caused  the  transmission  to  drop  at  low  energies. 
The  fiber  was  cut  into  two  parts  and  their  transmission  was  measured.  One  60  mm  long  part  also  has  low 
transmission  at  low  energies,  as  shown  in  figure  30.  It  was  suspected  that  the  channels  were  blocked  by 
glass  layer  of  thickness  t.  The  simulation,  with  wavi  ness  of  0.12  mrad  and  roughness  0.65  nm  obtained 
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in  single  fiber  simulation  previously,  and  adding  the  bending  with  radius  of  54  m,  fits  the  transmission 
spectrum  data  at  only  mid  and  high  energies.  Including  a  glass  layer  of  thickness  of  75  pm  into  the 
simulation,  the  simulated  transmission  fits  the  data  quite  well,  as  shown  in  figure  32.  Therefore,  the  drop 
in  transmission  at  low  energies  could  be  caused  by  glass  inclusions  in  the  channels.2 


Figure  33.  Sketch  of  fiber  displacements 
for  uniform  bending  measurements. 


Figure  34.  The  simulated  transmission  compared 
with  the  experimental  data  for  deliberately  bent 
fiber.  The  lines  are  simulated  data  with  waviness 
w=0.12  mrad,  and  roughness  z=  0.65  nm,  but 
with  different  bending  radius  R. 


In  the  straight  fiber  measurements,  as  shown  in  figure  8,  the  fiber  was  designed  to  be  straight  in  the 
experimental  setup.  However,  the  fiber  experienced  unintentional  bending  because  it  was  so  thin.  The 
bending  radius  was  determined  from  the  simulated  transmission,  fit  to  the  experimental  data.  In  the 
bending  experiment,  the  fiber  was  inserted  into  four  metal  screens,  as  shown  in  figure  33.  Initially,  the 
fiber  transmission  was  maximized  and  the  transmission  was  measured  and  analyzed  using  the  geometric 
simulation  as  for  the  straight  fiber.  The  bending  radius  found  from  the  simulation  was  called  initial  ra¬ 
dius  Ro  and  thus  produced  an  initial  fiber  displacement  from  an  ideal  straight  path.  Then  the  fiber  was 
bent  intentionally  by  moving  the  actuator  supporting  the  screen  as  in  figure  33.  The  nearest  screen  to  the 
input  end  of  the  fiber  was  fixed.  Each  screen  displacement  dmove  was  calculated  such  that  the  fiber  bent 
in  a  uniform  circle.  The  relation  between  the  bending  radius,  R,  and  the  screen  displacement,  dmovei,  is  R 
=  L2  /  2  dmovei,  where  Li  is  the  length  of  the  fiber  from  the  input  end  to  the  ith  screen.  In  the  initial  meas¬ 
urement,  the  fiber  had  been  bent  unintentionally  with  a  radius  Ro  which  resulted  unintentional  displace¬ 
ment  d0i ,  if  then  the  fiber  is  displaced  by  dmovei ,  the  total  displacement  of  the  fiber  at  the  ith  screen  is 
dj=  (d0i  +  dmovei  )■  A  fiber  with  length  of  140  mm  was  measured  in  the  bending  set  up.  The  initial  posi¬ 
tion,  the  fiber  assumed  to  be  straight,  was  found  from  simulation  to  have  an  unintentional  bend  with  ra¬ 
dius  35  m.  Deliberate  bends  of  R=9-23  m  were  measured.  The  simulation,  with  this  intentional  bending 
radius,  and  using  the  waviness  and  roughness  from  straight  fiber  simulation  (w=0.12  mrad  and  z=0.65 
nm),  fits  the  data  quite  well,  as  shown  in  figure  34.  The  good  agreement  between  the  computer  models 
and  the  single  fiber  data  allows  accurate  design  and  prediction  of  optic  behavior. 
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Figure  33.  Experimental  setup  for  transmission  measurement  of  lens. 


A  long  focal  length  multifiber  collimating  lens  with  a  length  of  140  mm,  input  size  of  29  mm  X  29 
mm  and  output  size  30  mm  x  30  mm  was  measured.  The  fibers  were  close-packed  at  the  input  end  of  the 
lens.  The  fiber  packing  fraction  at  the  output  end  of  the  lens  was  76.5%,  therefore  the  total  open  area  of 
the  lens,  including  the  individual  fiber  open  area  of  70%,  was  54%.  The  lens  was  tested  in  both  focus¬ 
sing  and  collimating  mode.  The  focal  length,  transmission  efficiency,  uniformity  and  exit  divergence 
were  measured  in  collimating  mode.  Scatter  rejection,  contrast  enhancement,  and  demagnification  was 
demonstrated  in  focussing  mode.  In  this  mode  the  detector  can  be  smaller  than  the  object  imaged. 

2.H.2.a  Collimating  mode 

The  experimental  setup  to  measure  focal  length,  transmission  efficiency,  and  uniformity  of  the  lens, 
is  shown  in  figure  35.  The  source  was  193  cm  from  the  input  of  the  lens,  the  detector  was  1  cm  from  the 
output  of  the  lens. 

2.H.2.a.l  Focal  length 

The  focal  length  of  the  lens  is  the  location  at  which  the  transmission  is  maximized.  Transmission  at 
different  distances  and  energies  are  shown  in  figures  36a  and  b.  The  measured  focal  length  of  the  lens 
was  193±1  cm.  An  image  of  the  x-ray  output  of  the  lens  with  this  source  distance  is  shown  in  figure  35. 

2.H.2.a.2  Transmission 

Transmission  as  a  function  of  energy  was  measured.  To  explain  the  experimental  results,  a  geometric 
simulation  for  single  fibers  was  extended  to  simulate  the  lens.  The  bending  radius  of  fiber  at  each  loca¬ 
tion  in  the  lens  was  calculated  using  a  simple  geometric  formula.  The  number  of  fibers  with  each  bend¬ 
ing  radius  was  determined.  Putting  together 
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Energy  (keV) 

Figure  34a.  Transmission  vs  energy  with 
difference  source-lens  distances.  The  trans¬ 
mission  is  maximized  at  1 92 ±1  cm. 


E=35keV 


Source-  Lens  Distance  (cm) 

Figure  34b.  Transmission  versus  source- 
lens  distance  at  various  energies.  For  all 
energies,  the  transmission  is  maximized  at 
a  192+1  cm. 


all  three  simulation  parameters,  waviness  and  roughness  from  the  single  fiber  simulation,  and  bending 
radius  calculated  at  each  fiber  location  in  the  lens,  the  simulated  transmission  spectrum  was  obtained. 
The  simulated  transmission  spectrum  for  the  whole  lens  is  the  average  transmission  spectra  of  all  the 
locations  normalized  with  the  total  fractional  open  area  of  the  lens.  The  simulated  transmission  spectrum 
compared  to  experimental  data  is  shown  in  figure  38.  The  simulated  transmission  fits  the  experimental 
data  better  at  higher  energies  than  lower  ones  since  the  channels  of  some  fibers  in  the  lens  were  blocked 
by  glass  inclusions  ,1’ 2  The  lens  has  good  transmission  up  to  40  keV,  and  therefore  seems  promising  for 
low  and  mid  range  energy  applications,  such  as  mammography,  astronomy  and  angiography. 

2.H.2.a.3  Uniformity 


Figure  37.  Photograph  of  the  output 
of  the  lens  with  the  source  at  the 
focal  distance. 
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Figure  38.  Simulated  transmission  com¬ 
pared  to  experimental  data  for  the  long 
focal  length  collimating  lens. 
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Figures  39  and  40  are  two-dimensional  output  scans  of  the  lens  at  20  keV  and  and  45  keV,  respectively. 
These  were  recorded  by  moving  a  5  mm  diameter  aperture  in  two  directions  perpendicular  to  the  lens 
axis, 


-10  -5  0  5  10  15 


Left-right 

Figure  39.  Output  intensity  from  two- 
dimensional  scan  of  lens  output  at  20 
keV  before  some  parts  of  the  fibers 
changed. 


Left-rigfrt(mm) 

Figure  40.  Output  intensity  from  two  di¬ 
mensional  scan  of  lens  output  at  45  keV 


with  source,  lens,  and  detector  fixed.  The  higher  the  energy,  the  lower  the  transmission  of  the  outer 
edges  of  the  lens.  This  is  expected,  because  of  smaller  critical  angle  at  higher  energies. 

2.H.2.a.4  Exit  Divergence 

Measurements  of  the  exit  divergence  of  the  collimating  lens  were  performed  at  17.5  keV  with  a  Mo 
source  and  Nal  detector  as  shown  in  figure  41.  The  total  divergence  was  measured  by  rocking  a  (100) 
silicon  wafer  in  the  beam,  using  the  (400)  Bragg  diffraction,  without  any  aperture  at  the  exit  of  the  lens. 
The  results  of  the  total  exit  divergence  are  shown  in  Figure  42.  Local  divergence  was  measured  with  a 
9.5  mm  diameter  aperture  at  upper  left  (ul),  upper  center  (uc),  upper  right  (ur),  center  left  (cl),  center  (c), 
center  right  (cr),  bottom  left  (bl),  bottom  center  (be),  and  bottom  right  (br)  on  the  exit  of  the  lens.  The 
results  of  the  local  divergence  are  shown  in  Table  III.  In  general,  the  divergence  is  low,  which  requires 
that  the  fibers  are  very  parallel  at  exit.  Some  areas  (uc  and  be)  had  poor  alignment,  and  the  fibers  were 
replaced.  The  average  local  divergence  is  in  agreement  with  the  value  of  the  total  divergence,  1 .2  mrad. 


Figure  41.  Experimental  setup  for  exit  divergence  of  lens 
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Local  divergence 

UL=1.2 

UC=1.5 

UR=1.0 

CL=0.9 

C=0.9 

CR=1.1 

BL=1.2 

BC=1.5 

BR=1.3 

Table  III.  Local  exit  divergence,  in 
mrad,  at  the  nine  aperture  locations. 


2.H.2.b  Focusing  Mode 

In  focusing  mode,  the  focusing  end  of  the 
lens  faced  the  detector,  thus  the  lens  acted  as  a 
collector  to  the  incoming  parallel  beam.  Testing 
of  this  lens  was  performed  at  the  Stray  Light  Fa¬ 
cility  (SLF)  at  Marshall  Space  Flight  Center 
(MSFC)  in  Huntsville,  Alabama  where  a  very  parallel  source  beam  was  available,  as  shown  in  figure  43. 
X  rays  from  a  copper  source  enter  a  vacuum  tube  100  m  in  the  length  and  exit  into  the  testing  area. 
Measurements  at  the  stray  light  facility  are  complicated  by  the  low  count  rate  of  the  photons  because  of 
the  huge  source-to-detector  distance.  Data  was  collected  over  a  10  keV  energy  window  15-25  keV,  to 
increase  the  count  rate. 


Crystal  Arge(degiee) 

Figure  42.  Total  divergence  rocking 
curve.  The  first  peak  is  diffraction  of  the 
Mo  Kai  and  the  second  the  Mo  Ka2 
characteristic  lines.  The  FWHM  of  the 
first  peak  was  1.2  mrad. 


2.H.2.b.l  Contrast  and  scatter  rejection  measurements 


Contrast  is  defined  as  the 
log  of  the  ratio  of  the  inten¬ 
sity  through  the  area  of  inter¬ 
est  to  the  intensity  of  the 
neighboring  background. 

C  =  In  ( Ip/Ibg)  (10) 

where  Ip  is  the  intensity 
through  object  of  interest,  and 
Ibg  is  the  neighboring  background  intensity.  Contrast  is  reduced  by  the  presence  of  scattered  x  rays 
which  emerge  at  all  angles  from  the  patient  and  form  a  background  fog.  Contrast  with  scatter 

Cs=ln{(Ip+Is)/(Ibg+Is)},  (11) 

where  Is  is  the  intensity  of  the  scatter. 

Scatter  fraction  is  defined  as 


x-ray 
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Figure  43.  Experimental  setup  for  contrast  and  scatter  measurement. 


SF=  ps/(Is+Ip)] 


(12) 
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Scatter  can  be  removed  from  the  beam  by  either  an  air  gap  or  a  polycapillary  optic.  The  introduction 
of  an  air  gap  causes  the  intensity  of  the  isotropic  scatter  radiation  to  drop  faster  than  the  forward- 
directed  primary  beam.  For  conventional  sources,  air  gaps  hurt  resolution  because  the  geometric  blur 
from  the  finite  source  increases  with  object  to  detector  distance.  Here,  the  parallel  beam  has  very  low 
geometric  blur.  Alternatively,  the  scatter  radiation,  which  emerges  from  the  phantom  at  high  angles,  can 
be  absorbed  by  the  glass  walls  of  the  polycapillary  optic,  as  shown  in  figure  44.  Only  the  forward- 
directed  primary  beam  is  carried  down  the  channels. 


Contrast  and  scatter  fraction  measure¬ 
ments  were  performed  in  focusing  mode  as 
shown  in  figure  41  with  a  detector-lens  dis¬ 
tance  of  20  cm.  The  scatter  fraction  was 
measured  by  placing  a  lead  block  in  the 
beam  and  measuring  the  scattered  fog  inten¬ 
sity  behind  the  block  due  to  a  2.5  cm  thick 
Lucite  phantom.  By  scanning  the  detector, 
the  scatter  intensity  behind  the  block  can  be 
measured,  as  shown  in  figures  45  and  46, 
with  parallel  air  gap  and  lens,  respectively. 
The  polycapillary  lens  is  even  better  than  the  parallel  beam  air  gap  at  removing  the  scatter.  To  measure 
contrast,  a  1  cm  thick  Lucite  phantom  with  holes  of  different  diameter  and  depth  was  attached  on  the 
output  of  the  vacuum  chamber.  Two  cases  were  measured.  In  the  first,  air  gap  case,  the  lens  was  re¬ 
moved  and  detector  was  placed  34  cm  from  the  phantom.  Because  the  air  gap  removes  most  of  the 
scatter,  the  contrast  is  good.  In  the  second  case,  the  lens  was  put  in  place  and  detector  inserted  20  cm 
from  the  lens.  Distance  between  the  phantom  and  the  lens  was  3  cm.  In  both  case,  the  output  was 
scanned  using  a  small  (3  mm  x  3  mm)  detector.  The  intensity  was  measured  for  60  seconds  at  each  lo¬ 
cation,  and  the  data  assembled  as  an  image,  as  shown  in  figures  47  and  48.  The  contrast  was  determined 
from  the  measured  intensity  behind  and  near  each  hole  and  is  shown  in  table  IV.  Contrast  was  also  cal¬ 
culated  from  the  theoretical  attenuation  in  the  plastic,  corrected  for  the  known  scatter  fractions.  In  each 


Figure  44.  Mechanism  photons  of  x  rays  through  the 
glass  wall 


Distance  from  lead  edge  (mm) 
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Distance  from  lead  edge  (mm) 


Figure  45.  Scatter  fraction  versus  distance 
from  lead  edge  with  parallel  beam  air  gap. 


Figure  46.  Scatter  fraction  versus  distance 
from  lead  edge  using  lens. 
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case  the  measured  value  agrees  quite  well  with  the  theoretical  calculation.  Because  the  scatter  fractions 
in  both  cases  are  low,  the  contrast  is  almost  the  same. 


Hole 

depth 

(mm) 

Experiment 

Theoretical 

Lens 

Parallel  beam 
air  gap 

Lens 

Parallel  beam 
air  gap 

4 

0.25±0.04 

0.24±0.03 

0.24 

0.23 

8 

0.48±0.04 

0.4110.03 

0.48 

0.47 

Table  IV.  Contrast  values,  measured  for  the  4  and  8  mm  deep 
holes  in  the  1  cm  thick  Lucite  phantom,  for  the  lens  and  air 
gap.  Theoretical  calculations  include  the  scatter  fraction  from 
figure  43  and  44. 
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Figure  47.  Image  created  from  a  Lucite  phantom 
with  different  diameter  and  depth  holes  with  a 
parallel  beam  air  gap. 


15 


20 


3 


.47 

.41 


059 


0  5  10 

XAxjs(imt 

Figure  48.  Image  is  created  from  a  Lucite  phantom 
with  the  lens. 


2.H.2.b.2  Demagnification 

The  image  in  figure  48  is  a  smaller  than  that  in  figure  47,  due  to  the  demagnification  produced  by  the 
focusing  of  the  lens.  The  amount  of  demagnification  can  be  calculated  from  the  change  in  center-to- 
center  distances  between  the  holes,  given  in  Table  V.  The  average  demagnification  is  0.87,  which  is  as 
expected  since  the  image  was  taken  approximately  nine  tenths  of  the  way  from  the  focal  point  of  the 
lens.  Demagnification  of  images  is  routinely  performed  in  digital  mammography  systems  by  converting 
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the  x  rays  to  visible  light  and  using  lenses  or  fiber  optics  to  decrease  the  size  of  the  image.  The  manipu¬ 
lation  of  the  x  ray  image  directly  may  allow  the  use  of  direct  digital  detectors,  which  do  not  require  a 
phosphor  for  conversion  to  visible  light,  and  therefore  can  have  higher  resolution. 

3.  Conclusion 

A  new,  more 
physical  model  for 
describing  profile 
defects  in  polycapillary 
x-ray  optics  at  high 
energies  for 
mammography  has 
been  developed.  The 
simulation  model  M2, 
is  physically  plausible 
and  works  quite  well  in 
explaining  the 
transmission  spectra.  It 

uses  a  clear  technique  for  determining  the  parameters,  which  makes  simulating  the  experimental  results 
and  designing  lenses  straightforward,  easy  and  timesaving.  According  to  the  simulations,  waviness  and 
bending  in  the  capillary  channel  profiles  can  be  particularly  harmful  for  high-energy  photons  because  of 
the  smaller  critical  angle  at  high  energies.  Roughness  is  not  important  in  the  simulation  for  borosilicate 
glass.  It  is  expected  that  increasing  the  rigidity  and  straightness  of  the  fiber  and  reducing  the  channel 
size  to  an  optimum  value  can  further  increase  the  transmission.  The  spectra  of  transmission  versus  en¬ 
ergy  have  been  found  to  be  a  very  sensitive  tool  in  capillary  quality  analysis,  and  transmission  versus 
source  scan  curve  measurements  can  also  provide  significant  information.  Lead  glass  shows  promise  in 
allowing  for  shorter,  higher  transmission  optics  which  will  easier  to  manufacture. 

The  collimating  lens  with  long  focal  distance  of  around  2  m  has  good  transmission,  uniformity  and 
low  exit  divergence.  The  good  agreement  between  the  computer  models  and  the  single  fiber  data  al¬ 
lowed  accurate  design  and  prediction  of  optic  behavior.  This  contributed  to  the  excellent  performance  of 
the  lens,  with  a  primary  transmission  of  almost  45  %  at  20  keV  and  nearly  zero  scatter  transmission, 
and  output  divergence  of  only  1 .2  mrad  at  17  keV.  Thus  the  lens  can  be  used  as  a  pre-patient  optic  to 
remove  the  variation  of  magnification  and  resolution  of  the  image,  and  as  a  post-patient  optic  to  reduce 
scatter  transmission  and  enhance  contrast  at  low  patient  dose.  Using  this  lens  in  focusing  mode  with  a 
parallel  beam  input,  demagnification  of  the  image  was  demonstrated.  The  lens  also  showed  equal  or 
better  scatter  rejection  and  contrast  compared  to  a  34  cm  air  gap. 

4.  Key  Research  Accomplishments 

•  A  new,  more  physical  model  for  describing  profile  defects  in  polycapillary  x-ray  optics  at  high  ener¬ 
gies  has  been  developed. 

•  A  geometric  simulation  with  two  fiber  defects,  waviness  and  bending  correction,  agrees  very  well  with 
experimental  data  for  borosilicate  fibers. 


Hole  center  to  cen¬ 
ter  distance  (mm) 

Lens 

Parallel 
beam  air  gap 

Demagnification 

ratio 

1-2 

6.9 

8.0 

0.87 

1-3 

9.0 

10.0 

0.89 

2-4 

9.6 

11.2 

0.86 

3-4 

6.9 

8.0 

0.87 

Table  V.  Change  in  distances  due  to  demagnification.  Average  demagnifl 
cation  is  0.87 ±0.02. 


Study  of  Capillary  X-ray  Optics  for  Digital  Mammography 


Wang,  Hui  and  Cari 


30 


•  The  same  geometry  simulation  with  four  fiber  defects,  waviness,  bending,  roughness  and  channel 
blockage,  agrees  very  well  with  experimental  data  for  leaded  glass  fibers. 

•  Measured  and  simulated  single  and  bent  borosilicated  and  lead  glass  fibers.  The  borosilicate  glass  had 
transmission  around  70  %  and  the  lead  glass  45  %  at  20  keV.  The  leaded  glass  has  better  scatter  rejetion 
than  borosilicate  glass.  Transmission  agree  well  with  simulations. 

•  Applied  simulation  from  single  fibers  to  a  multifiber  lens.  This  modeling  is  in  good  agreement  with 
measurements.  Therefore,  characteristics  from  single  fiber,  simulations  can  be  used  to  predict  the 
behavior  of  a  multifiber  lens. 

•  Measured  multifiber  lens  in  collimating  and  focusing  modes.  The  multifiber  lens  has  good  uniformity  ( 
non  uniformity  <  1.3  %) ,  high  transmission  ( 45  %),  and  low  divergence  (1.2  mrad).  This  makes  it  ef¬ 
fective  for  use  with  a  parallel  beam  input  as  a  post-phantom  demagnifying  optic.  The  demagnification  is 
0.87  and  agreed  with  prediction. 

•  The  multifiber  lens  demontrated  good  scatter  rejection,  which  occurs  because  of  the  tiny  critical  an¬ 
gles,  and  resulting  selective  reflection  of  on-axis  rays  by  the  glass  walls. 
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ABSTRACT 

Polycapillary  fibers  and  a  prototype  collector  for  high  energy  x  rays  with  a  2  m  focal  length  have  been 
fabricated  and  characterized.  Measurements  of  a  prototype  collector,  performed  in  collimating  mode,  show 
that  the  optic  has  high  transmission,  good  uniformity,  and  small  exit  divergence.  The  transmission  as  a 
function  of  energy  was  analyzed  using  an  extended  single  fiber  geometrical  optic  simulation  and  the  result 
shows  that  the  simulation  fits  the  data  fairly  well.  Scatter  transmission  and  contrast  enhancement  were 
measured  in  focussing  mode  using  a  parallel  beam  input. 


1.  INTRODUCTION 

Polycapillary  optics  can  be  used  to  improve  contrast  and  resolution  in  imaging.  X  rays  are  transported 
down  the  hollow  capillary  tubes  by  external  reflection  from  smooth  glass  walls.  A  single  polycapillary  fiber 
is  about  0.51  mm  in  diameter,  and  contains  hundreds  or  thousands  of  hollow  capillary  tubes,  10  pm  in 
diameter,  as  shown  in  figure  1 .  A  multifiber  lens  was  constructed  by  stringing  3039  flexible  singe  fibers 
through  metal  holding  grids.  The  fibers  at  one  end  of  the  lens  are  parallel,  and  at  the  other  end  are  pointed 
toward  a  common  focal  point,  2  m  away.  This  means  that  this  lens  can  be  used  to  collimate  or  focus  x  rays. 
In  collimating  mode,  the  lens  is  used  as  a  collimator  to  collect  x  rays  from  a  point  source  at  the  focus  and 
produce  a  nearly  parallel  beam.  In  focusing  mode,  the  lens  is  used  as  a  collector  for  parallel  a  beam  input 
and  produces  a  focused  beam.  The  focused  beam  can  be  used  to  demagnify  an  image.  This  is  of  potential 
benefit  for  direct  imaging  with  small  area  digital  x-ray  detectors. 


Figure  1.  SEM  of  a  single  fiber.  The 
diameter  of  the  channels  in  this  fiber  is 
about  50pm.  The  fibers  used  in  the 
experiment  have  open  area  of  70±  3  %, 
outer  diameter  of  0.51  mm  and  channel 
size  of  10  pm,  but  similar  structure. 


Figure  2.  3039  flexible  fibers  were  strung  through 
metal  grids  to  produce  the  multifiber  lens.  This  can  be 
used  as  either  a  collimating  or  focusing  lens.  One  end  is 
parallel  30  mm  x  30  mm.  At  the  other  end  (  29  mmx  29 
mm)  fibers  are  pointed  at  a  focal  point  2000  mm  away. 
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2.  BASIC  PRINCIPLE  OF  POLYCAPILLARY  OPTICS 

The  refractive  index  of  glass  for  x  rays  can  be  written  as 


n  =  J—  =  l-6-ip  (1) 

V  eo 

where  8  is  the  refractive  index  decrement  and  y  is  proportional  to  absorption  of  the  glass.  Because  the 
refractive  index  is  less  than  one,  x  rays  traveling  from  the  air  to  glass  will  be  totally  externally  reflected  at 
small  incidence  angles.  The  critical  angle  for  borosilicate  glass  is 


e,=V5s 


30mrad 

E(keV) 


where  0C  is  in  milliradians  and  E  is  the  photon  energy  in  keV.  For  example,  the  critical  angle  is  3.75  mrad 
for  8  keV  photons  and  0.375  mrad  for  80  keV  photons.  When  x-ray  photons  hit  the  glass  walls  of  the 
channels  at  angles  larger  than  critical  angle,  the  photons  will  be  absorbed.  Because  of  the  small  critical 
angle,  only  primary  x  rays  will  be  transmitted  through  the  capillary  channels  to  form  an  x-ray  image,  with 
near  total  scatter  rejection  of  Compton  scattered  x  rays. 


STRAIGHT  FIBER  MEASUREMENTS 


_100  cm.. 


Pinhole 


Sources  X-ray  beam 


Lead 

shielding 


Detector 


Figure  3.  Experimental  setup  for  single  fiber 
measurement 


-3-2-10123 

i  ,  ~  ,  .  Source  Position  (mm) 

The  experimental  setup  for  smgle  fiber  measurement  is 

shown  in  figure  3.  A  microfocus  source  with  a  tungsten  target  ^ure*.  ^e!?“rld  sourc* sca"  C"T  from 

and  source  spot  size  of  50  pm  was  used  in  the  experiment.  The  x-rays  leakage 

source  was  mounted  on  two  stages,  which  could  be  moved  in 

two  direction  transverse  to  the  beam,  and  also  enclosed  in  a  thick  lead  box  to  reduce  scattered  x-rays.  The 
distance  between  the  source  and  the  fiber  was  about  100  cm  to  ensure  incident  angles  less  than  critical 
angle  up  to  energy  80  keV.  A  200  pm  pinhole  was  attached  to  thick  lead  shielding.  The  size  of  the  pinhole 
is  smaller  than  the  diameter  of  the  fiber  which  is  about  0.5  mm  to  prevent  scattering  and  leakage,  to  reduce 
the  detector  dead  time  and  to  keep  the  distance  between  the  source  and  the  fiber  as  small  as  possible, 
because  of  critical  angle  for  high  energies.  The  fractional  open  area  of  the  fiber  is  70  ±  3  %.  The  fiber  was 
put  on  a  finely  machined  groove  in  an  aluminum  plate  mounted  on  two  stages,  which  could  be  moved  in 
two  directions  transverse  to  the  beam.  To  ensure  that  source,  pinhole,  fiber  and  detector  were  aligned  well, 
the  transmission  as  a  function  of  source  location  was  measured,  as  shown  in  figure  4  for  a  fiber  with  length 
of  140  mm.  The  figure  shows  that  there  was  not  any  x-ray  leakage  detected  and  the  pinhole  position  was  at 
the  center  of  the  fiber.  Then  the  transmission  of  the  fiber  type  used  to  construct  the  collimating  lens  were 
measured  in  the  10-80  keV  energy  range.  Several  fibers  of  different  lengths  were  measured,  three  fibers 
from  which  are  presented  in  figure  5. 
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3.1  Simulation  Analysis 
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Figure  5.  Transmission  versus  energy  for  single 
fibers  with  three  different  lengths.  The  lines  are 
simulated  data  with  waviness,  w=  0.12  mrad,  and 
roughness,  z=  0.65  nm,  but  with  different  bending 
radius,  R140=210  m,  R125=105  m,  R60=48  m. 
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Figure  6.  Simulated  transmission  compared 
with  the  experimental  data  for  a  defected 
fiber.  Best  fit  is  simulation  with  waviness, 
w=0.12  mrad,  roughness,  z=0.65  nm, 
bending  radius,  R=54  m,  and  glass  layer 
thickness,  t=75  pm. 


The  experimental  data  then  was  analyzed  using  a  geometrical  optic  simulation  to  understand  the 
performance  of  the  fibers.  The  geometric  algorithm  is  a  two  dimensional  approximation.  The  geometric 
simulation  includes  three  fiber  defects:  bending,  waviness  and  roughness1.  Bending  is  the  deviation  of  the 
channel  from  an  ideal  straight  path,  approximated  as  uniform  bending  with  a  radius  R.  Waviness  is  local 
tilting  of  channel  walls,  parameterized  by  a  standard  deviation  of  tilt  angles,  w.  Surface  roughness  is  short 
wavelength  surface  variation,  parameterized  as  an  r.m.s.  displacement  z.  The  best  fit  simulation 
parameters  were  determined  in  the  manner  of  reference  1 .  Using  waviness  from  the  source  scan  curve  at  30 
keV,  roughness  from  source  scan  curve  at  10  keV  and  adding  a  bending  radius  of  210  m,  a  fit  was  achieved 
across  the  whole  spectrum.  The  simulated  source  scan  curves  with  these  three  best  fit  simulation 
parameters  for  a  140  mm  long  fiber,  compared  with  the  experimental  data,  is  shown  in  figure  5.  Then  the 
fitting  processes  were  applied  to  two  fibers  with  length  of  125  mm  and  60  mm.  The  simulated  transmission 
spectra  for  these  three  fibers  compared  with  experimental  data  are  also  shown  in  figure  5.  Each  simulation 
fits  the  experimental  data  well.  Since  the  fibers  are  thin,  a  shorter  fiber  bends  more  and  therefore  transmits 
somewhat  more  poorly  at  high  energies.  Ali  three  lengths  of  fibers  shown  in  figure  5  have  good 
transmission  up  to  energy  50  keV. 


3.2  Defected  Fiber 


One  of  the  fibers  measured  contained  defects  which  caused  the  transmission  to  drop  at  low  energies. 

The  fiber  was  cut  into  two  parts  and  their  transmission  was  measured.  One  part  with  length  of  60  mm  also 
has  low  transmission  at  low  energies,  as  shown  in  figure  6.  It  was  suspected  that  the  channels  were  blocked 
by  glass  layer  of  thickness  t.  The  simulation,  with  waviness  of  0.12  mrad  and  roughness  0.65  nm  obtained 
in  single  fiber  simulation  previously,  and  adding  the  bending  with  radius  of  54  m,  fits  the  transmission 
spectrum  data  at  only  mid  and  high  energies.  Including  a  glass  layer  of  thickness  of  75  pm  into  the 
simulation,  the  simulated  transmission  fits  the  data  quite  well,  as  shown  in  figure  6.  Therefore,  the  drop  in 
transmission  at  low  energies  could  be  caused  glass  inclusion  in  the  channels.2 
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Figure  7.  Sketch  of  fiber  displacements  for 
uniform  bending  measurements. 


Energy(keV) 

Figure  8.  The  simulated  transmission  compared 
with  the  experimental  data  for  deliberately  bend 
fiber.  The  lines  are  simulated  data  with  waviness, 
w=0.12  mrad,  and  roughness,  z =  0.65  nm,  but  with 
different  bending  radius,  R. 


In  the  straight  fiber  measurements,  as  shown  in  figure  3,  the  fiber  was  designed  to  be  straight  in  the 
experimental  setup.  However,  the  fiber  experienced  unintentional  bend  because  it  was  so  thin.  The  bending 
radius  was  determined  from  the  simulated  transmission,  fit  to  the  experimental  data.  In  the  bending 
experiment,  the  fiber  was  inserted  into  four  metal  screens,  as  shown  in  figure  7.  Initially,  the  fiber 
transmission  was  maximized  and  the  transmission  was  measured  and  analyzed  using  the  geometric 
simulation  as  for  the  straight  fiber.  The  bending  radius  found  from  the  simulation  was  called  initial  radius 
Ro  and  thus  produced  an  initial  fiber  displacement  from  an  ideal  straight  path.  Then  the  fiber  was  bent 
intentionally  by  moving  the  actuator  supporting  the  screen  as  in  figure  7.  The  nearest  screen  to  the  input 
end  of  the  fiber  was  fixed.  Each  screen  displacement  dmove  was  calculated  so  that  the  fiber  bent  in  a 
uniform  circle.  The  relation  between  the  bending  radius,  R,  and  the  screen  displacement,  di?  is  R  =  L*  /  2  di? 
where  Li  is  the  length  of  the  fiber  from  the  input  end  to  the  im  screen.  Since  in  the  initial  measurement  the 
fiber  had  been  bent  unintentionally  with  a  radius  R0,  and  then  displaced  dmove ,  the  total  displacement  of  the 
fiber  at  the  i*  screen  is  d*=  (doi  +  dmove ).  A  fiber  with  length  of  140  mm  was  measured  in  the  bending  setup. 
The  initial  position,  the  fiber  assumed  to  be  straight,  was  found  from  simulation  to  be  unintentional  bend 
with  radius  of  23  m.  The  simulation,  with  the  bending  radius  calculated,  as  above,  and  using  the  waviness 
and  roughness  from  straight  fiber  simulation  (w=0.12  mrad  and  z=0.65  nm),  fits  the  data  quite  well,  as 
shown  in  figure  8.  The  good  agreement  between  the  computer  models  and  the  single  fiber  data  allows 
accurate  design  and  prediction  of  optic  behavior. 


4.  MEASUREMENT  AND  ANALYSIS  OF  A  LONG  FOCAL  LENGTH  MULTIFIBER 

COLLIMATING  LENS 

A  long  focal  length  multifiber  collimating  lens  with  a  length  of  140  mm,  input  size  of  29  mm  X  29  mm 
and  output  size  30  mm  x  30  mm  was  measured.  The  fibers  were  close-packed  at  the  input  end  of  the  lens. 
The  fiber  packing  fraction  at  the  output  end  of  the  lens  was  76.5%,  therefore  the  total  open  area  of  the  lens 
including  the  individual  fiber  open  area  of  70%  was  54%.  The  lens  was  tested  in  both  focussing  and 
collimating  mode.  The  focal  length,  transmission  efficiency,  uniformity  and  exit  divergence  were  measured 
in  collimating  mode.  Scatter  rejection,  contrast  enhancement,  and  demagnification  were  demonstrated  in 
focussing  mode.  In  this  mode  the  detector  can  be  smaller  than  the  object  imaged. 
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4.1  Experimental  Setup 

4.1.1  Collimating  mode 

The  experimental  setup  to  measure  focal  length,  transmission  efficiency,  and  uniformity  of  the  lens,  is 
shown  in  figure  9.  The  source  was  193  cm  from  the  input  of  the  lens,  the  detector  was  1  cm  from  the  output 
of  the  lens. 


4.1. 1.1  Focal  length 

The  focal  length  of  the  lens  is  the  location  at  which  the  transmission  is  maximized.  Transmission  at 
different  distances  is  shown  in  figure  10.  The  measured  focal  length  of  the  lens  was  192±1  cm.  An  image  of 
the  x-ray  output  of  the  lens  with  this  source  distance  is  shown  in  figure  1 1 . 


4.1.1.2  Transmission 

The  transmission  as  a  function  of  energy  was 
measured.  To  explain  the  experimental  result,  a 
geometric  simulation  for  single  fibers  was 
extended  to  simulate  the  lens. 

The  bending  radius  of  fibers  at  each  location  in 
the  lens  was  calculated  using  a  simple  geometric 
formula.  The  number  of  fibers  with  each  bending 
radius  was  determined.  Putting  together  all  three 
simulation  parameters,  waviness  and  roughness 
from  the  single  fiber  simulation,  and  bending 
radius  calculated  at  each  fiber  location  in  the 
lens,  the  simulated  transmission  spectrum  was 
obtained.  The  simulated  transmission  spectrum 
for  the  whole  lens  is  the  average  transmission 
spectra  of  all  the  locations  normalized  with  the 
total  fractional  open  area  of  the  lens.  The 
simulated  transmission  spectrum  compared  to 
experimental  data  is  shown  in  figure  12.  The 
simulated  transmission  fits  the  experimental  data 
better  at  higher  energies  than  lower  ones  since 

the  channels  of  some  fibers  in  the  lens  were  blocked  by  glass  inclusion  2  The  lens  has  good  transmission 
up  to  40  keV,  and  therefore  seems  promising  for  low  and  mid  range  energy  applications,  such  as 
mammography,  astronomy  and  angiography. 


Figure  10.  Transmission  vs  energy  with  difference 
source-lens  distances.  The  transmission  is  maximized  at 
193  cm. 


4.1. 1.3  Uniformity 

Figures  13  and  14  are  two-dimensional  output  scans  of  the  lens  at  45  keV  and  and  60  keV,  respectively. 
These  were  recorded  by  moving  a  5  mm  diameter  aperture  in  two  directions  perpendicular  to  the  lens  axis, 
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Figure  11.  Photograph  of  the  output  of  the 
lens  with  the  source  at  the  focal  distance. 


Energy  (keV) 


Figure  12.  Simulated  transmission  compared  to 
experimental  data  for  the  long  focal  length 
collimating  lens. 


with  source,  lens,  and  detector  fixed.  The  higher  the  energy,  the  lower  the  transmission  of  the  outer  edges 
of  the  lens.  This  is  expected,  because  of  smaller  critical  angle  at  higher  energies. 
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Figure  13.  Output  intensity  from  two 
dimensional  scan  of  lens  output  at  45  keV 


Figure  14.  Output  intensity  from  two 
dimensional  scan  at  60  keV 


4.1. 1.4  Exit  Divergence 

Measurements  of  the  exit  divergence  of  the  collimating  lens  were  performed  at  17.5  keV  with  a  Mo 
source  and  Nal  detector  as  shown  in  figure  15.  The  total  divergence  was  measured  by  rocking  a  (100) 
silicon  wafer  in  the  beam,  using  the  (400)  Bragg  diffraction,  without  any  aperture  at  the  exit  of  the  lens. 
The  results  of  the  total  exit  divergence  are  shown  in  Figure  16.  Local  divergence  was  measured  with  a  9.5 
mm  diameter  aperture  at  upper  left  (ul),  upper  center  (uc),  upper  right  (ur),  center  left  (cl),  center  (c),  center 
right  (cr),  bottom  left  (bl),  bottom  center  (be),  and  bottom  right  (br)  on  the  exit  of  the  lens.  The  results  of 
the  local  divergence  are  shown  in  table  I.  In  general,  the  divergence  is  low,  which  requires  that  the  fibers 
are  very  parallel  at  exit.  Some  areas  (uc  and  be)  had  poor  alignment,  and  the  fibers  were  replaced.  The 
average  local  divergence  is  in  agreement  with  the  value  of  the  total  divergence,  1.2  mrad. 
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Figure  15.  Experimental  setup  for  exit  divergence  of  lens 


Crystal  Angle  (degree) 


Local  divergence 

UL=1.2 

UC=1.5 

UR=1.0 

CL=0.9 

C=0.9 

CR=1.1 

BL-1.2 

BC=1.5 

BR=1.3 

Table  I.  Local  exit  divergence,  in  mrad  at 
the  nine  aperture  locations. 


Figure  16.  Total  divergence,  rocking  curve.  The  4.2  Focusing  Mode 

first  peak  is  diffraction  of  the  Mo  Koti  and  the 

second  the  Mo  Ka2  characteristic  lines.  The  In  focusing  mode,  the  focusing  end  of  the  lens 

FWHM  of  the  first  peak  was  1.2  mrad.  faced  the  detector,  thus  the  lens  acted  as  a 

collector  to  the  incoming  parallel  beam.  Testing  of 
this  lens  was  performed  at  the  Stray  Light  Facility  (SLF)  at  Marshall  Space  Flight  Center  (MSFC)  in 
Huntsville,  Alabama  where  a  very  parallel  source  beam  was  available.  X  rays  from  a  copper  source  enter  a 
vacuum  tube  100  m  in  the  length  and  exit  into  the  testing  area.  Measurements  at  the  stray  light  facility  are 
complicated  by  the  low  count  rate  of  the  photons  because  of  the  huge  source-to-detector  distance.  Data  was 
collected  over  a  10  keV  energy  window  15-25  keV,  to  increase  the  count  rate. 


4.2.1  Contrast  and  scatter  rejection  measurements 


x-ray  beams 


Vacuum  chamber 


< 


100  m 


Contrast  is  defined  as 
the  log  of  the  ratio  of  the 
intensity  through  the  area 
of  interest  to  the  intensity 
of  the  neighboring 
background. 

C  =  log  ( Ip/Ibg)  (3) 


Figure  17.  Experimental  setup  for  contrast  and  scatter  measurement. 


where  Ip  is  the  intensity 
through  object  of  interest, 
and  Ibg  is  the  neighboring 
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background  intensity.  Contrast  is  reduced  by  the  presence  of  scattered  x  rays  which  emerge  at  all  angles 
from  the  patient  and  form  a  background  fog.  Contrast  with  scatter  is 


Cs=ln(Ip+Is)/(Ip+Is)5  (4) 

where  Is  is  the  intensity  of  the  scatter. 

Scatter  fraction  is  defined  as 

SF=  [IS/(IS+IP)]  (5) 

Scatter  can  be  removed  from  the  beam  by  either  an  air  gap  or  a  polycapillary  optic.  The  introduction  of  an 
air  gap  causes  the  intensity  of  the  isotropic  scatter  radiation  to  drop  faster  than  the  forward-directed  primary 
beam.  For  conventional  sources,  air  gaps  hurt  resolution  because  the  geometric  blur  from  the  finite  source 
increases  with  object  to  detector  distance.  Here,  the  parallel  beam  has  very  low  geometric  blur. 
Alternatively,  the  scatter  radiation,  which  emerges  from  the  phantom  at  high  angles,  can  be  absorbed  by  the 
glass  walls  of  the  polycapillary  optic,  as  shown  in  figure  18.  Only  the  forward-directed  primary  beam  is 
carried  down  the  channels. 

Contrast  and  scatter  fraction 
measurements  were  performed  in 
focusing  mode  as  shown  in  figure 
17  with  a  detector-lens  distance  of 
20  cm.  The  scatter  fraction  was 
measured  by  placing  a  lead  block 
in  the  beam  and  measuring  the 
scattered  fog  intensity  behind  the 
block  due  to  a  2.5  cm  thick  Lucite 
phantom.  By  scanning  the 
detector,  the  scatter  intensity  behind  the  block  can  be  measured  as  shown  in  figure  19  and  20.  The 
polycapillary  lens  is  even  better  than  parallel  beam  air  gap  at  removing  the  scatter.  To  measure  contrast, 
a  1  cm  thick  Lucite  phantom  with  holes  of  different  diameter  and  depth  was  attached  on  the  output  of  the 
vacuum  chamber.  Two  cases  were  measured.  In  the  first,  air  gap  case,  the  lens  was  removed  and  detector 
was  placed  34  cm  from  the  phantom.  Because  the  air  gap  removes  most  of  the  scatter,  the  contrast  is  good. 
In  the  second  case,  the  lens  was  put  in  place  and  detector  inserted  20  cm  from  the  lens.  Distance  between 
the  phantom  and  the  lens  was  3  cm.  The  output  was  scanned  using  a  small  (3  mm  x  3  mm)  detector.  The 
intensity  was  measured  for  60  seconds  at  each  location,  and  the  data  assembled  as  an  image,  as  shown  in 
figures  21  and  22.  The  contrast  was  determined  from  the  measured  intensity  behind  and  near  each  hole  and 
is  shown  in  table  II.  Contrast  was  also  calculated  from  the  theoretical  attenuation  in  the  plastic,  corrected 
for  the  known  scatter  fractions.  In  each  case  the  measured  value  agrees  quite  well  with  the  theoretical 
calculation.  Because  the  scatter  fractions  in  both  cases  are  low,  the  contrast  is  almost  the  same. 


0.18 
0.16 
0.14 
O  0.12 

2  010 
t  0.08 

1  0.06 

o 

w  0.04 
0.02 
0.00 

0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0 
Distance  from  lead  edge  (mm) 


Distance  from  lead  edge  (mm) 


Figure  18.  Mechanism  photons  of  x  rays  through  the  glass  wall. 


Figure  19.  Scatter  fraction  versus  distance  from  Figure  20.  Scatter  fraction  versus  distance  from 
lead  edge  with  parallel  beam  air  gap.  lead  edge  using  lens. 
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4.2.2  Demagnification 


The  image  in  figure  22  is  a 
smaller  than  that  in  figure  21,  due 
to  the  demagnification  produced  by 
the  focusing  of  the  lens.  The 
amount  of  demagnification  can  be 
calculated  from  the  change  in 
center-to-center  distances  between 
the  holes,  given  in  Table  III.  The 
average  demagnification  is  0.87, 
which  is  as  expected  since  the  image  was  taken  approximately  nine  tenths  of  the  way  from  the  focal  point 
of  the  lens.  Demagnification  of  images  is  routinely  performed  in  digital  mammography  systems  by 
converting  the  x-  ray  to  visible  light  and  using  lenses  or  fiber  optics  to  decrease  the  size  of  the  image.  The 
manipulation  of  the  x-ray  image  directly  may  allow  the  use  of  direct  digital  detectors,  which  do  not  require 
a  phosphor  and  conversion  to  visible  light,  and  therefore  can  have  higher  resolution. 


Hole 

depth 

(mm) 

Experiment 

Theoretical 

Lens 

Parallel  beam 
air  gap 

Lens 

Parallel  beam 
air  gap 

4 

0.25+0.04 

0.24+0.03 

0.24 

0.23 

8 

0.48+0.04 

0.41+0.03 

0.48 

0.47 

Table  II.  Contrast  values,  measured  for  the  4  and  8  mm  deep  holes  in 
the  1  cm  thick  Lucite  phantom,  for  the  lens  and  air  gap.  Theoretical 
calculations  include  the  scatter  fraction  from  figure  19  and  20. 


XAxis  (mm) 

Figure  21.  Image  created  from  a  Lucite 
phantom  with  different  diameter  and  depth 
holes  with  a  parallel  beam  air  gap. 
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Figure  22.  Image  is  created  from  a  Lucite 
phantom  with  the  lens. 


■ijjjmwmigg 

Lens 

Parallel  beam 
air  gap 

Demagnification 

ratio 

1-2 

0.87 

1-3 

0.89 

2-4 

11.2 

3-4 

HEH 1 

mammm 

Table  III.  Change  in  distances  due  to  demagnification.  Average 
demagnification  is  0.87±0.02. 


5.  CONCLUSION 

The  good  agreement  between  the  computer  models  and  the  single  fiber  data  allowed  accurate  design  and 
prediction  of  optic  behavior.  This  contributed  to  the  excellent  performance  of  the  lens,  with  a  transmission 
of  >  38  %  at  10-40  keV,  and  output  divergence  of  only  1.2  mrad  at  17  keV.  Using  this  lens  in  focusing 
mode  with  a  parallel  beam  input,  demagnification  of  the  image  was  demonstrated.  The  lens  also  showed 
equal  or  better  scatter  rejection  and  contrast  compared  to  a  34  cm  air  gap. 


a*  * 
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